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Abstract. 

Anyons are quasiparticles in two-dimensional systems that show statistical 
properties very distinct from those of bosons or fermions. While their isolated 
observation has not yet been achieved, here we perform a quantum simulation of anyons 
on the toric code model. By encoding the model in the multi-partite entangled state of 
polarized photons, we are able to demonstrate various manipulations of anyonic states 
and, in particular, their characteristic fractional statistics. 
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1. Introduction 

In three spatial dimensions, only two types of statistical behaviors have been observed 
dividing particles into two groups: bosons and fermions. If one is restricted to two- 
dimensional systems the situation changes. There, anyons [1] can appear, which exhibit 
fractional statistics that ranges continuously from bosonic to fermionic. Anyons are 
responsible for the fractional quantum Hall effect [2] and it has been demonstrated 
that they could be realized as quasiparticles in highly entangled many-body systems. 
Formally, the properties of anyons are described by two-dimensional topological 
quantum field theories [3] that dictate their trivial dynamical, but complex statistical 
behavior. In general, it is expected that such topological quantum field theories come 
into effect at low energies of highly correlated many-body systems, such as quantum 
liquids p[]. However, the observation of anyonic features requires high population in the 
system's ground state, high purity samples and, above all, the ability to separate the 
anyonic effects from the dynamical background of the strongly interacting system. In 
spite of significant experimental progress [5] the fractional statistics of anyons has not 
yet been conclusively confirmed due to the complexity of these systems [6]. 

To demonstrate characteristic anyonic features we employ here a different strategy. 
We simulate the toric code model [7] in a quantum system without the continuous 
presence of interactions [Hj. The toric code is a two-dimensional topological lattice 
model, where anyons are spatially well localized. Recently, several proposals showed how 
the toric code could possibly be implemented on extended lattices of qubits enabling 
one to employ anyonic states for, e.g., protected quantum computation [9]. While 
such a quantum computational scheme requires large systems, here we show that a 
number of anyonic properties can be demonstrated by considering four-partite entangled 
Greenberger-Horne-Zeilinger (GHZ) states [10] . These states are dynamically encoded 
in the polarization of photons, rather than produced via cooling [8J. This provides the 
additional advantage of being subject to negligible decoherence due to the weak coupling 
of photons to the environment. Polarized photons as qubits are a well understood and 
controllable quantum optical system that already allowed to observe the four-qubit GHZ 
state [TT] . Here we employ this system to simulate anyonic behavior. 

2. Toric code model 

2.1. Hamiltonian and ground state 

The anyonic model under consideration is based on the toric code proposed by Kitaev [7J 
and is described in detail in tutorial introductions [12] • This system can be defined on a 
two-dimensional square lattice with interacting qubits placed at its vertices. To facilitate 
the description of the interactions let us split the lattice into two alternating types of 
plaquettes labelled by p and s, as in figure [Q The defining Hamiltonian is 

IV \ „z _z _z „_z \ " _a; _x /i \ 

H = ~ 2s a p,l (T P ,2 (T p,Z (T pA ~ 2s <J s,l (T s,2 (T s,3 (T sA > I 1 ) 
p s 
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where the summations run over the corresponding plaquettes and the indices 1, ...,4 of 
the Pauli operators, a z and a x , enumerate the vertices of each plaquette in a counter- 
clockwise fashion. Each of the s- or p-plaquette interaction terms commute with the 
Hamiltonian as well as with each other. Thus, the model is exactly solvable and its 
ground state is explicitly given by [13] 

\0 = n ^(1 + <:<,<X ;!(")...()) , (2) 

with o- z \0) = |0). The state |£) represents the anyonic vacuum state and it is unique for 
systems with open boundary conditions. 
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Figure 1. The toric code lattice with qubits at the vertices of a square lattice 
with two types of plaquettes. Light blue and red plaquettes correspond to the s- 
and p-plaquettes, respectively. Qubit rotations enable manipulations of anyons on 
neighboring plaquettes. (a) Application of a z on a single qubit yields two e-type 
anyons placed at the neighboring s plaquettes, where the string passes trough the 
rotated qubits. Similarly, m anyons are created on p plaquettes by a x rotations, (b) 
Two a x rotations create two pairs of m-type anyons. If one anyon from each pair 
is positioned on the same plaquette then they annihilate, thereby connecting their 
strings, (c) When a part of a string forms a loop around unpopulated plaquettes, the 
loop cancels (dashed), (d) Here we restrict to one s plaquette and four neighboring p 
plaquettes. Anyon e is produced by a a z on qubit 1 (blue), |e) = of This system 
can support the circulation of an m around an e anyon. 



2.2. Anyonic quasiparticles 

Starting from this ground state one can excite pairs of anyons connected by a string 
on the lattice using single qubit operations. More specifically, by applying a z on some 
qubit of the lattice a pair of so called e-type anyons is created on the two neighboring 
s plaquettes [figure [Qa)] and the system is described by the state |e) = o" z |£)- An 
m pair of anyons lives on the p plaquettes and is obtained by a a x operation. The 
combination of both creates the composite quasiparticle e with |e) = a z a x \£) = io y \£). 
Two equal Pauli rotations applied on qubits of the same plaquette create two anyons 
on this plaquette. The fusion rules (lxl = exe = mxm = exe = l,exm = e, 
1 x e = e, etc., where 1 is the vacuum state [12]) describe the outcome from combining 
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two anyons. In the above example, if two anyons are created on the same plaquette then 
they annihilate. This operation also glues two single strings of the same type together 
to form a new string, again with a pair of anyons at its ends [Figure [11(b)]. If the string 
forms a loop, the anyons at its end annihilate each other, thus removing the anyonic 
excitation. In case that only a part of the string forms a loop, the string gets truncated 
[Figure (He)]. For non-compact finite systems, such as the one we consider here, a string 
may end up at the boundary describing a single anyon at its free endpoint. 

2.3. Anyonic statistics 

Anyonic statistics is revealed as a non-trivial phase factor acquired by the wave function 
of the lattice system after braiding anyons, i.e., after moving an m anyon around an e 
anyon [Figure [0(d)] or vice versa. Consider the initial state l^ini) = of|£) — l e )- If an 
anyon of type m is assumed to be at a neighboring p plaquette it can be moved around 
e along the path generated by successive applications of a x rotations on the four qubits 
of the s plaquette. The final state is 

= crX^X'l^ini) = -af{afa^afat |£» = -|*w> ■ (3) 

Such a minimal loop, which vanishes the moment it is closed, is analogous to 
the application of the respective interaction term C s = l* 7 ® 2 a s 3 a s 4 ( or = 
a p,i a p,2 a p,2, a p,A) °f ^ ne Hamiltonian. This operator has eigenvalue +1 for all plaquettes 
of the ground state whereas it signals an excitation, e.g., | \J/ ini ) , with eigenvalue 

— 1, when applied to the plaquette where an anyon resides. However, ([3]) is much more 
general, as the actual path of the loop is irrelevant. It is the topological phase factor of 

— 1, which reveals the presence of the enclosed anyon. Alternatively, we can interpret 
(J3j) as a description of twisting e, the combination of an e and an m-type anyon, by 
2tt. The phase factor of —1 thereby reveals its 47r-symmetry, characteristic for half 
spin, fermionic particles [H] . Note that the e (m) anyons exhibit bosonic statistics with 
respect to themselves [TJ. 

3. Experimental implementation 

3.1. Minimal instance of the toric code model 

As this two-dimensional system is well suited for demonstrating characteristic anyonic 
features, the question arises how big the lattice has to be. It turns out that four qubits 
of a single s plaquette represent the minimal unit of the toric code model [TS] . In 
this case, the four neighboring p plaquettes are represented only by their corresponding 
links. Hence, while the presence of an e anyon is detected by the plaquette operator 
C s , the presence of m anyons is detected by the eigenvalue of the four C' p = o z v i a^^ 
operators that correspond to the adjacent links of the p plaquettes. Consequently, the 
product determining the ground state [Eqn. (j2J)] reduces to only one factor for the single 
s-plaquette resulting in a four qubit GHZ state, | £) = (| 0000) + | This is an 
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eigenstate of the relevant C s and C operators with eigenvalue +1. An eigenvalue — 1 
for any of these operators denotes the presence of the corresponding anyon. 




nii<<<<mx<<<< 
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Figure 2. (a) The experimental set-up for the observation of the anyonic vacuum state 
of the toric code model, the demonstration of anyonic features and the verification of 
anyonic statistics. The photons are created in a 2 mm thick (3 Barium Borate (BBO) 
crystal, which is pumped by femtosecond ultraviolet (UV) pulses. Walk-off effects are 
compensated for by a HWP and a 1 mm thick BBO crystal. For the observation of the 
anyonic states, we are interested in the second order emission of the SPDC [16], initially 
emitted in the two spatial modes a and b: cx |2H) Q |2V) 6 + |2V) a |2H) h + |HV) a |HV) h . 
The photons are further processed by a HWP in mode b, which transforms |H) into 
|+) and |V) into |-) [with |±) = (|H) ± \V))/V2}. The state |£) is observed behind 
the PBS and the BSs under the condition of detecting a photon in each of the four 
modes 1, 2, 3, 4 (average count rate: 2.8 min -1 ). The local operations (LO) for the 
demonstration of anyonic features are implemented by operation specific HWPs and 
quarter- wave plates (QWPs). The polarization state of each photon is analyzed (PA) 
with a HWP and a QWP in front of a PBS. (b) Measurement outcomes of |£) in the 
basis of for each qubit. 



3.2. State implementation 

In our simulation the qubits supporting the anyonic states are encoded in the 
polarization of single photons propagating in well-defined spatial modes. This means 
that for the anyonic vacuum and for all obtained final states our system is of the form 
|GHZ*) = (|H 1 H 2 H 3 H 4 ) + e i ^|V 1 V 2 V 3 V4))/v / 2. The indices (omitted in the following) 
label the spatial mode of each photon, i.e., they represent the position of the qubit as 
in Figure LUd), and H (V) denote linear horizontal (vertical) polarization, representing 
a logical (1). To obtain such four-photon entangled states, the second order emission 
of a non-collinear type-II spontaneous parametric down conversion process (SPDC) 
[17] . yielding four photons in two spatial modes, is overlapped on a polarizing beam 
splitter (PBS) and afterwards symmetrically split up into four spatial modes by two 
polarization independent beam splitters (BS), see Figure [2]^a) [IB] . Prior to the second 
order interference at the PBS the polarization of two of the photons is rotated by a 
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half- wave plate (HWP). Under the condition of detecting one photon in each spatial 
mode we observe the desired states. 

As we are interested in the vacuum |£) = |GHZ°) and the anyonic state |e) = 
<7 i\0 = I GHZ 7 "), the successful simulation and demonstration of anyonic features relies 
on a careful distinction and characterization of these two orthogonal GHZ states. 
This comprises the confirmation of genuine four-partite entanglement and a method 
to reveal the phase 0. For this purpose, state analysis is performed by measuring 
the correlation, c z , in the of basis for each qubit and the correlation function c xy {^) = 
0i (7) <g) 02 (7) $03(7)® 04 (7) with 0^(7) = [(cos 7) of + (sin 7) erf] [19]. This measurement 
is ideally suited to prove the coherent superposition of the terms |HHHH) and |VVVV) 
and to determine their relative phase 0. It is obtained from a measurement of all qubits 
along the same direction in the x-y plane of the Bloch sphere. For a GHZ^ state its 
expectation value shows an oscillation, (0^(7)) = Vcos (47 + 0), with a period of one 
quarter of 2tt. This is a unique signature of four qubit GHZ entanglement. The visibility 
of this oscillation V is 1 for pure GHZ states. For a general state it equals twice the 
element Phhhh.vvvv — Tr(|HHHH)(VVVV|p) of the corresponding density matrix p. 

To analyze the performance of the set-up, first, the observation of the state 
|£) = I GHZ ) is confirmed. Figure Wljo) shows all possible measurement outcomes 
in the 0i02 <T 3 cr 4 basis, where the expected populations of |HHHH) and |VVVV) are 
clearly predominant with probabilities of Phhhh = (41.2 ± 3.4)% and Pvvvv = 
(39.6 ± 2.7)%, respectively, close to the expected Phhhh = -Pvvvv = 50%. We find 
a correlation of (c z ) = 0.908 ± 0.047. The population of other terms is caused by noise 
arising from higher-order emissions of the down conversion and a remaining degree 
of distinguishability of photons at the PBS. The dependence of (c xy {^)) for the state 
|£) is displayed in Figure [3](a) from which we infer a visibility V of (68.3 ± 1.1)%. 
This value is obtained from a weighted least squares fit to a Fourier decomposition of 
( c xy{l)) considering only even components up to order 4: {c xy {pf)) = J2t=o a k cos(2A; • 
7) + bksin(2k ■ 7). Only these components can originate from physical states. From 
the phase of the correlation function we deduce the phase of the GHZ state to be 
= (0.02 ± 0.01)77, close to the expected value of 0. 

The visibility V together with the populations Phhhh and Pvvvv allow further to 
determine the fidelity F = (GHZ|p|GHZ) = (V + Phhhh + P V vvv)/2. We obtain a 
value of (74.5 ± 2.2)% (for measurement data see also Table [[]). The fidelity is not 
only important for an estimation of state quality but can also be applied to verify 
genuine four-partite entanglement, an essential element in the presented toric code 
model. With a proper entanglement witness a fidelity greater than 50% is sufficient [20J. 
The experimentally observed fidelity is clearly above this bound, i.e., it is of high enough 
quality to allow for the demonstration of the anyonic properties. 
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3.3. Anyonic manipulations and detection of anyonic statistics 

We start by analyzing the state change under the creation (and annihilation) of anyons, 
thereby demonstrating the characteristic fusion rules. Applying o\ in mode 1 creates 
an e-type anyon on the s plaquette resulting in a GHZ 71 " state, Figure [3](b). This is 
clearly proven by the phase <fi = 1.027T of the correlation function. A further application 
of aj on any other mode, say j = 3, changes the state of the plaquette according to 
e x e = 1, and the anyon is moved away from the s plaquette under investigation. 
These two a z rotations represent a string connecting two anyons, which traverses the 
particular plaquette without influencing its state, as confirmed by the observation of 
the initial GHZ state [(f) = O.OI71", Figure E^c)]. A further crj rotation on one of the 
remaining vertices, e.g., j = 4, creates an e occupation on the plaquette. Thus, the 
observation of GHZ 71 [0 = 1.037T, Figure EJ^d)] demonstrates the non-trivial version of 
the 1 x e = e fusion rule. Alternatively, this can be seen as demonstrating the invariance 
of an anyonic state, when a string traverses the plaquette. 
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Figure 3. Pictorial representations of various evolutions of anyons, together with the 
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(a) The vacuum state, (b) Generating an e-type anyonic state |e) = IGHZ^) by 
a a z operator, (c) A second a z removes the anyon from the plaquette and elongates 
the string resulting in the GHZ state, (d) A third a z creates another pair of anyons, 
connected with the string already traversing the plaquette. (e) Moving an m anyon 
around an empty s plaquette or (f) first creating an e anyon on the plaquette (i), 
performing the loop (ii), and removing the e anyon from the plaquette again (iii) give 
identical final states up to a global phase factor, (g) The global phase, it, is revealed 
by superimposing the two above evolutions, thereby clearly proving anyonic statistics, 
(h) An alternative path for the anyon e that gives the same fractional statistical phase, 
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To detect the major feature of anyons, their non-trivial statistical phase acquired 
when moving one anyon around the other, we employed an interference measurement, 
which makes their overall phase factor visible [see ([3]) and Figure ED(d)]. Let us study 
the two evolutions separately. The first is to create a pair of two m anyons, move 
one of them around the empty s plaquette and then annihilate them. This evolution is 
equivalent to having a string ending up at a p plaquette with an m anyon at its endpoint 
and circulating that anyon around the s plaquette. The obtained loop operator C s , 
obviously, does not change the vacuum state and, thus, results in GHZ [<f> = — 0.027T, 
Figure [3](e)]. Alternatively, we (i) create an e anyon on the plaquette, (ii) encompass 
it with the loop of an m anyon, and (iii) remove it again [Figure 09(f)] • The whole 
evolution is, in analogy to Q, described by <j%{<j\ <rf cxf o"f)crf|£) = — 1£) = — |GHZ°). 
The correlation function determines faithfully the angle = — 0.0l7r [Figure 19(f)], but 
is blind to the characteristic overall phase factor. This, finally, can be observed by 
(i) first generating the superposition e -i7r/,4 (|£) + i\e))/y2 of having an anyon e on the 
plaquette or not by the unitary operation (cxf) _1//2 ; (ii) moving the m anyon around this 
superposition gives e -Mr / 4 (|£) — z|e))/y2, the superposition of the above evolutions; and 
(iii) the application of the inverse rotation (of) 1//2 makes the phase difference visible 
resulting in —i\e) = — ^GHZ 71 "). The observed value of = 1.007T [Figure E9(g)] clearly 
proves the phase acquired by the braiding and therefore the anyonic statistics [21] . 
To reinforce the observation of the obtained fractional phase, we perform an additional 
interference experiment with an alternative path. Now, the superposition of an e anyonic 
state and the ground state is generated in the s-plaquette by rotating one qubit and it 
is removed by rotating another qubit, as seen in Figure [3](h). If we move an m anyon 
around the plaquette in between the two steps, we observe again the characteristic phase 
shift of -1. 

4. Conclusions 

Our results show that we can create, manipulate and detect the toric code states by 
encoding them in a simple physical system. In our quantum simulation we demonstrated 
several of the fusion rules and we provided supporting evidence for the fractional 
statistics of the toric code anyons. Multiqubit toric code states are known to be useful 
for novel types of quantum error correction [22] . Extending our experimental results 
presented here to larger, scalable quantum systems [231 [21] will enable the application 
of the toric code for quantum information processing in the future [25] . 

Note added: With the completion of this work we became aware of a similar 
experimental implementation, closely resembling |8], that was simultaneously performed 
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a z x a z 2 \i) (0.02 ± 0.01) -vr 


(77.4 ±2.8)% 


a\al\i) (0.01 ± 0.01) • tt 


(77.3 ± 2.5)% 


a z o-l\£) (-0.01 ± 0.01) -tt 


(76.3 ± 2.6)% 


String passing through the s plaquette populated with an anyon 


<7f<7f|e) (1.02 ± 0.01) -vr 


(75.2 ± 2.4)% 


cr|(7| e) (1.03 ± 0.01) • tt 


(76.7 ±2.7)% 


<jf(7fe) (1.02 ± 0.02) -tt 


(74.5 ±3.1)% 


Loop around an unpopulated s plaquette 




C s \0 (-0.02 ± 0.02) ■ vr 


(76.6 ± 3.4)% 


Loop around a populated s plaquette followed by annihilation of the anyon 


{<Jl)C s {<Jl)\i) (-0.01 ± 0.01) • tt 


(76.8 ± 2.8)% 


Interference procedure to reveal anyonic statistics 




K) 1/2 CsK)~ 1/2 IO (1.03 ± 0.01) • tt 


(76.1 ±3.0)% 


(aiy^Cslal)- 1 / 2 ^) (0.99 ± 0.01) • it 


(73.5 ± 3.6)% 


K^Of)-!/^) (1.01 ± 0.01) • TT 


(75.2 ±3.0)% 


K) 1/2 ^K)~ 1/2 I0 (1.00 ± o.oi) • tt 


(75.8 ± 2.5)% 


to) 1 ' 2 c a te)- 1 ' 2 \t) (l.oo ± o.oi) -vr 


(78.7 ±3.2)% 



Table 1. Experimental measurement data of the anyonic states. 
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